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Using Fermat's least optical path principle the famif^-^j: ray-trajectories through a special but 
common type of a gradient refractive index lens, n (r) ^-^o + AnR/r, is solved analytically. The 
solution, i.e. the ray-equation r (</>), is shown to be closely related to the famous Rutherford scat- 
tering and therefore termed photonic Rutherford scattering. It is shown that not only do these 
classical limits correspond, but also the wave-mechanical^ictures coincide: The time-independent 
Schrodingier equation and the inhomogeneous Helmholz elat ion permit the same mapping between 
massive particle scattering and diffracted optical scalar waves. Scattering of narrow wave-packets 
finally recovers the classical trajectories. The analysis suggests that photothermal single particle 
microscopy infact measures photonic Rutherford scattermg^in specific limits. 



I. INTRODUCTION 



Almost exactly 100 years ago in the year 1911, Ernest 
Rutherford changed our picture of the atom by his fa- 
mous theory on the scattering of positively charged a- 
particles^. In Rutherford scattering positively charged 
Helium nuclei are deflected by a Coulomb potential orig- 
inating from positive nuclei of gold atoms as originally 
shown by Rutherford, Geiger and Mardsen 2 -". This work 
has been a milestone in the discovery of the structure of 
the atom, revealing that most of the mass of an atom 
is concentrated in a tiny nucleus. Thus, Rutherford 
scattering is considered in each atomic physics lecture, 
treated in a classical framework to provide the char- 
acteristic angular distribution of scattered a particles. 
While a classical showpiece illustrating Rutherford scat- 
tering may be obtained from a paraboloidal hard wall- 
potentiaP, a direct display of the continuous trajectory 
or measuring a single deflection instead of the total cross- 
section remains difficult. The classical theoretical pre- 
dictions by E. Rutherford were later revisited to account 
for the detailed structure of the atom. While small im- 
pact parameters could be used to systematically probe 
the core-potentiaP, large impact parameters needed to 
additionally account for the electronic shielding of the 
core Coulomb potential. Somewhat unexpectedly^, the 
intricate 8 quantum-mechanical spin-less treatment of the 
Coulomb 1 / r-potential predicted for all energies the same 
scattering-cross section as the classical theor^. Here we 
present the photonic analog of Rutherford scattering. It 
is given in the geometrical optics approximation (GOA) 
by the deflection of rays (the classical limit) or, in wave 
optics, as the diffraction of waves by a 1/r-refractive in- 
dex profile. This profile is provided by a heat point- 
source in a homogeneous medium. Such a point source 
may be a light- absorbing nano-particle 10 embedded in 
some medium w hich a re used in photothermal single par- 
ticle microscop ^ 11 * 12 *. Experimental dem onstrations of 
the effect can be achieved (see Section IV). 



prese&ded and in Section [TIT] the analytical solution de- 
rived^In Section [IV] analogies of the found ray-optics 
solutiCTTare explored with respect to the classical non- 
relativtstic and relativistic Rutherford scattering prob- 
lem ^hout radiation reaction. In Section [V] the wave- 
mecka*iical pictures are explored. Here, the correspon- 
dence^etween QM Coulomb scattering and the scalar 
opticaAn field in the 1/r-profile inhomogeneous refractive 
indexffield is revealed. Thereafter the correspondences 
to tr^-dassical pictures are established. Both an optical 
Fresno-diffraction and a QM wave-packet formalism are 
used t£ achieve the necessary departure from the plane- 
wav^Hpit. Finally, the found solutions are applied to 
phol^^ermal microscopy and compared to previous ex- 
periiMe^its. 
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The paper is structures as follows: In Section [H] the 
ray-optics treatment of the 1 / r-refractive index profile is 
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II. CLASSICAL LIMIT: FERMATS' PRINCIPLE 

Obtainable through a variational principle with fixed 
path end-points which unifies Maupertius' (mechanics) 
and Fermat's (optics) variational principle, the following 
differential equation suitable for massive particles and 
light may be obtained ^ 3 1 14 1 23 !: 



ds 2 



V ( ^n 4 (r)v 2 (r) 





dr 


)• 





n(r) v(r) , (1) 



with r being a vector on and s a stepping parameter along 
the path. The difference in treating light or massive par- 
ticles consists in the proper choice of the velocity v(r). 
In the latter case one may take n = 1, such that Eq. 
([I]) reduces to Newton's first law, Eq. (J2|, and thus also 
classical dynamics with the choice of the stepping pa- 
rameter ds = d£, by setting v 2 /2 = E/m — V/m, i.e. the 
specific difference of total and potential energy per unit 
maseP. Eq. (l) may even be used to describe relativistic 
gravitational mechanics in a static space-ti me m etric by 
its corresponding non-unit refractive indexLUiLLl] To de- 
scribe the paths of rays of light, Eq. is to be supple- 
mented by setting v = c/n, where c is the vacuum speed 
of light. This case will correspond to Fermat's princi- 
ple of the least optical path and allows the calculation of 
light trajectories through a spatially inhomogeneous re- 
fractive index field n (r). This picture provides a classical 
particle picture of light propagation and corresponds to 
the zero- wavelength limit of wave-optics^. The result, 
Eq. (J3|, is the "F=ma" -optics developed by Evans et al. 
and explored by many others^"^ 



mechanics : 
optics : 



d^r 



W(r), 



dr 
di 
dr 



= «(r), (2) 
n(r). (3) 



While the solution of positive energies to the Newton's 
equation of motion, Eq. (pi, on a 1/r-potential is known 
as Rutherford scattering (see Section IV), we will now 
seek the physically achievable analogon in the optical 
domain. Consider a heat source that generates a tem- 
perature profile T (r) = T + AT (r) with 



AT(r) = P abs /(4^r), 



(4) 



which, according to Fouriers law, decays with the inverse 
distance r from the object to To at infinite distance (P a bs 
and k are the absorbed power and the medium heat con- 
ductivity, respectively). This temperature profile results 
in the linear regime in the refractive index profile Eq. ([5| 
that takes up the inverse distance dependence with the 
thermo-refractive coefficient dn/dT as a proportionality 
factor, 



n (r) = n + ^ AT (r) = n + An 



R 



(5) 



where no = n (To) is the unperturbed real- valued refrac- 
tive index, R the radius of the heat-source and An = 



AT (R) dn/dT a real- valued refractive index contrast. 
This is valid as long as the thermal conductivity of the 
finite-size heat source is larger than the mediums' con- 
ductivity. As we will demonstrate, the problem of find- 
ing the ray-trajectories fulfilling Eq. ([3| is equivalent 
to the scattering by an unshielded Coulomb potential, 
i.e. Rutherford scattering. A similar but rather artificial 
type of refractive index field, n 2 (r) = const. + 2/c/r, has 
been shown to yield all types of Kepler-orbits for light 
in that mediunP^S Also, effective refractive indices 
have been shown to mimic the path of light in gravi- 
tational fields as pr edicted by Einsteins theory of general 
relativit ^ 13 | 14 | 18 | 2Q -"^. For the weak gravitational field 
limit of the Schwarzschild metric n (r) = l + 2GMc~ 2 r~ 1 
describes the null geodesies of light. 



III. EXACT SOLUTION 

Since, by symmetry, the trajectories will be con- 
fined to a plane (see Fig. [T]), we use cylindrical co- 
ordinates (r, (j)) where the acceleration takes the form 
r" = r(r" -rcj)' 2 ) + (r<j>" + 2r'<j>') and the gradient 
reads Vn = f d r n + 6r~ 1 don = n _1 Vn 2 /2. The prime 
denotes differentiation with respect to the stepping pa- 
rameter s. Fermats' least optical path principle Eq. (|3| 
then gives two equations, Eq. (|6| for the radial coordinate 
and Eq. ^ for the angular coordinate: 

f : r" - r<j)' 2 -n AnR-^ 



-An 2 R z - (6) 

/ - r" J 

attractive/repulsive attractive 

: r<t>" + 2r>'= (7) 

The above set of coupled differential equations is equiv- 
alent to the perturbed Kepler problem with its precess- 
ing orbit solution^. Equation (JT]) yields the conserved 



optical angular momentum L z 



r>2 AJ 



, i.e. L' 



0. 



Now, the formula of Bouguer 16 allows to express this 
quantity at any point along the trajectory as L z = 
rsin((/>) |dr/ds|, such that with Eq. (|3| we find at in- 
finite distance L z = briQ. The parameter b > 0, the so 
called impact parameter, is the distance of the approach- 
ing parallel ray to the optical axis (see Fig. [TJt). The dif- 
ferential Eq. ^ is the analogue to the mechanical radial 









b- 










FIG. 1. Annot ated sketch of an exemplary ray trajectory 
(red) r (0), Eq. (12), through the refractive index fields n (r) 
with (a) An > and (b) An < in Eq. (5}. 
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force equation and shows an inverse radius squared inter- 
action, which is either attractive or repulsive depending 
on the sign of An, and a perturbation by a inverse ra- 
dius cubed term (underlined in the following). To solve 
it for r (0), a change of differentials is needed. Applying 
d/ds = (j) f d/d(j) = L z r~ 2 d/d(j) twice, and introducing the 



inverse radius variable u 
relation 



1/r one finds the following 



2„,2 



ds 2 



Liu 



;_d_l 

d(j) \~ d(j) u 



2, ,2 1 



-L„u 



la' 



which transforms Eq. (pi) into 



L 2 u 2 



d 2 u 



Liu 3 



-n 



AnRu 2 -An 2 R 2 u 3 



We now introduce the variable 

no 



AnR' 



(8) 



(9) 



(10) 



which is a measure for the inverse strength of the heat 
induced refractive index gradient and encodes the polar- 
ity of the interaction in such a way that a positive sign of 
£ corresponds to repulsion. Equation J9| then becomes, 
after rearranging and collecting of the terms linear in u, 



"u 

la 



(ii) 



If the refractive index in the medium is homogeneous, 
i.e. £ = oo, the harmonic oscillator differential equation 
with unit angular frequency emerges and the correct so- 
lution fulfilling the boundary conditions is u = r _1 = 
6 -1 sin ((/>). In cartesian coordinates y = rsin(</>) this is 
a straight line parallel to the optical axis at a distance 6, 
which of course is the unperturbed ray, see dashed line in 
Fig. [TJl. If the perturbation is nonzero, and requiring for 
the moment that > 1, Eq. (11) has the form of the 



familiar harmonic oscillator differential equation plus a 
constant, u" + c\u = — C2 with positive c\. It is solved by 

- l] with the yet to be de- 
Equation (11) is therefore 



u = ^ [e cos {y/c[ ((j) - O )) 
termined constants e and </>o 
solved by 



P 



e cos 



(7[0-<AoD- 



(12) 



where eccentricity is allowed to be either positive or neg- 
ative, and with the parameters 



P =[ve-i\/i 

7 2 =1-6-2^-2 



(13) 



Mathematically, the orbits described by Eq. (12) repre- 



sent perturbed hyperbolic trajectories with the particle 
being the exterior (£ > 0) or interior (£ < 0) focuJ^, 
see Fig. [lji,b. More exactly, they are epispirals, a special 
case of so-called Cotes's spirals. Such orbits may show 
peculiar behavior, such as multiply revolving trajectories 



for £ < 0, when the perturbation-parameter 7 approaches 
zero (see Fig. [2]) and were already discussed by the grand- 
son of Charles Darvin, C. G. Darwin, in the context of 
relativist ic Rutherford scattering of electrons in 191*P^. 
Also, somewhat later in 1916, Sommerfeld in his relativis- 
ts corrections to the Hydrogen spectra encountered the 
bound form of such orbits for the electrorP^U. To ob- 
tain the eccentricity e we reconsider the particular choice 
of the stepping parameter in Eq. J3|, and write again in 
cylindrical coordinates: 



— » 



+ rV = n(r) z 



(14) 



The radius of closest approach is obtained by setting 
r' = 0, and yields r m = b + £ _1 . Again, angular mo- 
mentum conservation (j)' = L z r~ 2 was used. Comparison 
of this expression to the corresponding minimum radius 
as described by Eq. (12), r m = pj (e — 1) at the angle of 
closest approach </> = 0o, yields the eccentricity e = b£. 
Setting the denominator of Eq. Jl2| ) to zero yields the 
extreme angles 6^ = ±|7 _1 | arccos (1/e) + 0o- Requir- 
ing that the ray approaches parallel to the optical axis 
from negative infinity, see Fig. |TJd, i.e. = 7r, will ori- 
ent the solution Eq. ( 12 ) according to the imposed initial 



conditions. We then find the angle of closest approach: 



= 7r — I7 1 | arccos (e 



(15) 



The parameters in Eqs. (13) and ([19]) together with Eq. 



(12) now fully determine the ray-trajectory. The scatter- 
ing angle = 6^, i.e. the deflection angle of an incoming 
horizontal ray, may be expressed as 6 = 20o — tt- We 
finally note that the differential Eq. Jill) can also be ob- 
tained from Binet's orbit equatio n 1 24 1 ^ 1 with the correct 
identification of the force terms as given by Eq. J3|. 

The previous treatment relied on the assumption, 
which is however valid in practical situations, that > 
1. If the impact parameter gets very small, 7 would be- 
come imaginary. This situation is solemnly due to the 
presence of the attractive inverse cubic interaction term 
which dominates the inverse squared term at small dis- 
tances, see Eq. J6|. Instead of Eq. (11), we must then 
solve the following differential Eq. : 



d(t>- 



(b~ 2 c 2 



1 



(16) 



It has the form u" — c\u = —c^ with positive c\ and is 
solved by u = ^ [ecosh (y/c\ (4> — 4>o)) + l] , where we 
have chosen the nyperbolic cosine for now and will con- 
sider the general solution hereafter. Therefore, 



r r (cj)) = 



Pr 



. COsh (j r [(j) - 0o,r]) + 1 ' 



(17) 



with the positive perturbation parameter j r > 1 de- 



termined by % 



-7 and p r 



[l — b 2 ^ 2 ] /£ = —p. The only admittable solution for an 
approach from infinity is for an eccentricity to be within 
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— 1 < e r < 0. In this situation r m>r = b + £ _1 > for 
£ > (and only for the repulsive case) is achieved at 
4> = 0o an d indeed yields e r = — 6£ in the desired range. 
Setting the denominator of Eq. ( fl7| ) to zero one finds the 
extreme angles 6^ r = ±|7~ 1 |arccosh (fr -1 £ -1 ) +</>o such 
that again we correctly orient the solution with the choice 
^oo r = 71 an d thereby </>o, r = tt — |7~ 1 |arccosh (fr -1 £ -1 ). 
Here, too, the deflection angle is r = 2</>o — tt ~ 
7T + 26f ln(6f/2) + (9 (6 2 £ 2 ) and its limit is (9 H> tt as 
6£ — 0, which corresponds to a perfect retroreflection for 
a head-on impact of a ray onto the lens. Both {0o,r? #r} 
are smooth continuations of {(/>o, 0} found earlier. Infact, 
allowing the cosine to have a complex argument with 
7 r = in Eq. (12), the same solution is obtained. 



We now seek the general solution of Eq. (16): 



r a (<t>) 



Pa 



e ,i exp (j a <j>) + e a> 2 exp (-j a <j>) + 1 ' 



(18) 



This ansatz now allows the incoming ray to have the 
correct distance at infinity, e.g. lim^-^ sin ((/>) r s ((/>) = 6, 
and gives the set of two two generalized eccentricities: 



e ,i 
e a ,2 



2&7a 



(19) 



Solution (18) works for both the attractive case, hence 
the subscript a , and the repulsive case. In the former 
case the solution is a true mixture of the hyperbolic sine 
and cosine which describes trajectories approaching from 
infinity and falling within a finite time into the coordi- 
nate origin. It does so without a closest distance r m and 
coming from the bt; < — 1 case the rays revolve evermore 
vigorously around the origin. Both phenomena continue 
the limiting behavior of Eq. (12) where the closest ap- 



proach distance goes to zero and the scattering angle 6 
diverges to infinity as \b£\ — » 1, see Fig. (J2|. In the case 
of repulsive interaction the solution given above guises 
the solution involving only the hyperbolic cosine found 
earlier, i.e. Eq. (17). 



Similar to certain cases of relativistic point-particle 
Kepler mechanics^, these two special solutions involv- 
ing the hypergeometric functions corresponds to special 
types of Cotes spirals. 



IV. PHOTONIC RUTHERFORD SCATTERING 

As the refractive index change itself is typically small 
for most materials (| An| « 10 -3 ), and since b > R for the 
incoming rays, the product ^> 1 is a large number. 
This allows allows us to approximate Eq. ( 12 ) by 



r{4>) 



V / PFTTcos(0-0 o )±l 



(20) 



where ± is the sign opposite of £, which now shows 
complete equivalence to the classical (non-relativistic) 



Rutherford scattering solution of Eq. Q on the potential 
V (r) = CV" 1 (attractive: C < 0, repulsive: C > 0), 



^RF (0) 



2Eb 2 /\C\ 



± 1 



(21) 



where the notion is such that attractive interaction is 
represented by the upper and repulsive interaction by 
the lower sign, respectively. The scattering parameters 
are 



E = mvl/2 

C = qiq 2 / (47re ) 

e 2 = 4E 2 b 2 C~ 2 -hi, e > 

(j)o = 7r ± arccos (1/e) 



(22) 



and describe the total energy E and mass m of the scat- 
tered particle, and e 2 the squared eccentricity of the or- 
bit, C the Coulomb force constant for the two charges 
#i 5 2 that describes the mechanical force F (r) = —Cr~ 2 r. 
The scatterer is assumed to be fixed here, i.e. has an in- 
finite mass as compared to the scattered particle. The 
angular momentum of the particle relative to the scat- 
terer at the origin is L = mvob, while its specific angular 
momentum (twice the areal velocity) is L z = L/m. The 
deflection of photons by a weak gradient index lens gen- 
erated by a heated point-like absorber, described by Eq. 
(20), is thus the complete photonic analogon of Ruther- 



ford scattering of a particles on a single nucleus, Eq. (21 ). 
V — » — n (r) 2 / '2 +^q/2 w Uq £ -1 r -1 can therefore be iden- 
tified as the photonic analogon of the potential energy 
decaying to zero at infinite distance, E — >• Uq/2 being 
the total energy and C —> —u^AnR is the equivalent of 
the Coulomb force constant, as can be inferred from Eq. 
(|6|. The form of Eq. ([3| also requires the mass to be 
set to unity ra = 1 in optics. Hence, all further equa- 
tions, e.g. the differential scattering cross section (^) 
unravelling the famous sin -4 (#/2) dependence, or the to- 
tal cross-section a>e of scattering by an angle larger than 
some angle can be obtained using these equivalences 
and the substitution 2E/C — >• £. 

The observation that the total energ y is positive re- 
quires a few comments. Typicall y 14 * 16 * it is stated that 



Eq. 



¥ 2 



0, corresponds to the equation for 



the total energy analogon, comprised of a kinetic energy 
term \ and a potential energy term — |n 2 , and thus 
the total energy in the optical case is equivalent to the 
mechanical scenario at zero energy E — 0. However, due 
to the inclusion of the additional constant shift (+71q/2) 
of the potential energy scale in V, necessitated by includ- 
ing no m Eq. (|J), we find that here the mechanical zero- 
energy scenario does not represent the optical problem at 
hand (cf. footnote 15 of the Evans et al. "F=ma" -optics 
paper 16 ). Indeed, the parabolic unit-eccentricity orbits of 
zero-energy scattering is not the found (approximate) so- 
lution for the ray-trajectory. Here, ^|r'| 2 + V — tiq/2 = 
and E can be identified with the first two terms yielding 
E = Uq/2 to be taken as the mechanical total energy 
analogon. Thus, only the unbound (hyperbolic) trajec- 
tories from classical mechanics are attainable for no ^ 0. 
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The discrepancy by a factor of 2 between the expres- 
sion for the distance of closest approach tq and the exact 
value from Eq. (14), r m (0) = stems from the fact 
that for b — >• the validity of the approximation bt; ^> 1 
and thus Eq. (20) breaks down. For a repulsive potential 
Eq. ( fl7| ) then passes the point of closest approach. The 
same argument explains the difference between r m i n (b) 
and r m (b). For the repulsive case (An < 0), Eq. (|6| 
shows an additional attractive inverse radius-cubed in- 
teraction resulting in a closer approach. Solving Eq. ( [l4| ) 
without such a term yields the photonic r m i n (fr)-value 
listed in Table [IJ The exact trajectories will thus pene- 
trate the classical Rutherford shadow region given by the 
paraboloid r s = 4£ -1 / [1 — cos ((/>)]. For large b^ ^> 1 the 
two expressions coincide. 



100 




FIG. 2. Absolute scattering/deflection angle \0\ (left axis) and 
the normalized distance of closest approach r m in (b) /b vs. im- 
pact parameter b for fixed interaction strength Black 
dashed-solid lines: Rutherford scattering, red lines: exact 
solution (orange dashed: attractive). Clearly visible is the 
effect of the additional attractive perturbative force allowing 
closer approaches and weaker deflections for the repulsive case 
(£ > 0) and stronger deflection in the attractive case (£ < 0). 
For large bt; ^> 1 both results converge. 

The before-men tioned similarity to relativistic motion 
in a 1/r-potentia l 25 ' 28 ' 29 ' (cf. paragraph §39 of reP^), 
which was used by Arnold Sommerfeld to give the 
fine-splitting of the Hydrogen line-spectrum 26 , may be 
brought to a correspondence with the photonic problem 
here using the previous substitutions complemented by 
the additional rule c — ^ no and with E = mc 2 — > n 2 , 
now replacing the total energy including the rest-mass. 
However, in optics there is no such distinction between 
relativistic and non-relativistic treatments of light, such 
that it shall not be implied that optics corresponds to me- 
chanics in its special relativistic form. Still, the analogy is 
between non-relativistic mechanics and optics as embod- 
ied in Eq. Furthermore, these relativistic Rutherford 
scattering solutions should not be taken as necessarily 
being more accurate since also here radiation reaction 
as embodied in the Lorentz-Abraham-Dirac equation are 
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TABLE I. Correspondence table showcasing different expres- 
sions in photonic and Rutherford/Coulomb scattering. 



not considered (see Huschilt et al! 31 * 32 * or Aguiar et alP^ 
and references therein). 

Rutherford scattering is generally considered for mul- 
tiple scattering on many nuclei with random impact pa- 
rameters. Thus, the measurable cross-section delivers 
the same results for attractive or repulsive Coulomb in- 
teractions. While the latter results are obtained by a 
classical and QM wave- mechanics, similar prediction for 
diffraction on multiple refractive index profiles have been 
made^. The photonic equivalent can, however, be tested 
easily on a single scattering center, allowing to access 
even the sign of the interaction, i.e. the sign of dn/dT, 
with the help of simple photodetector. For this, a macro- 
scopic experiment with a metal sphere embedded in a 
transparent resin as shown in Fig. [4] can be setup. Upon 
heating of the central sphere by a high-power laser one 
can measure the deflection of paraxial thin laser-beams 
according to cot (0/2) = b£. Also, quite naturally the 
photonic Rutherford scattering can be seen by the un- 
aided eye directly. Viewing an object through such a 
medium containing a heat point-source will cause the 
viewed object to appear warped according to the extrap- 
olated path as seen in Fig. [3] As noted before, a refractive 
index profile o f n (r) = 1 + 2GMc~ 2 r~ 1 describes gravi- 
tational lensinj^ 14 ' 18 1 20 ^ . Therefore, the observed dis- 
tortion nicely model for instance the famous Einstein ring 
phenomenon if a material with An > is used (some 
types of glasses such as N-PK51 have this property, cf. 
the TIE19 data sheet "Temperature Coefficient of the Re- 
fractive Index" of the manufacturer Schott). A computer 
program interactively visualizing these effects is publicly 
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available on the authors' webpage. 

Photothermal single particle microscopy also pro- 
vide measurements on single photonic Rutherford 
scatterers^l. A simple formalism starting from the ray- 
optics results presented here have been used to provide a 
semi-quantitative minimal model for photothermal lens- 
ing microscopy of heatable metallic nanoparticles 35 . 
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FIG. 3. Visual effect of the photonic scatterer (laser- heated 
metal sphere in transparent resin). An image photographed 
through the medium is warped. The deflection of rays gives 
the illusion of a crunching of the image. Warped image a) 
and the initial image b) were mirrored along the horizontal 
due to the layering turbidity from the manufacturing process. 
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V. WAVE MECHANICAL RUTHERFORD 
SCATTERING 

Similar to scattered alpha particles, also photons obey 
the wave-particle duality. While they interact only very 
weakly among each other, their interaction with mat- 
ter is described in its strength by the dielectric function 
e. The dielectric function thus defines a "photonic po- 
tential" manipulating the propagation of photons or op- 
tical scalar fields in the simplest form of wave optics. 
Similarly, Quantum Mechanics is the wave-description of 
matter. While the equivalence of the wave-optics treat- 
ment of the Coulomb scattering problem with the clas- 
sical has been shown for the plane- wave case, we will 
here demonstrate the correspondence also in optics and 
with finite beams (of either particles, or rays). Concepts 
from scalar wave optics have been applied to Quantum 
problems ever since and show the close re lation of both. 
As seen in nuclear scattering experiment J 36 * 37 *, molecu- 
lar interferrometry datsP^, or atomic aperture diffraction 
experiments 39 , interference effects for instance may con- 
veniently be described by Fresnel diffraction. Further 
mappings have been found between paraxial wave-optics 
and the Schrodinger equation in two dimensions'"^. We 
will now show the equivalence between our recent scalar 
wave optics treatment of the diffraction by the refrac- 
tive index profile n (r), Eq. (J5|, and the QM problem of 
scattering on a bare Coulomb potential. 



A. Plane wave scattering / diffraction 



The Schrodinger equation (SE) for scattering on a 
Coulomb potential V (r) = Cr~ l reads — |^;V 2 \I/c + 
-*c = E^c = hmvZSfrc. The wave vector k of the 



>c = u^c = 2 rm; o ^co- 
incident particle-wave is defined through the de Broglie 
relation Tik = ttivq. The time-independent SE may then 
be rewritten as 



2uk 
r 



(23) 



FIG. 4. Macroscopic Experiment on single Rutherford- like 
photonic scatterer (black-body sphere with a small hole). For 
a typical polymer medium, k ~ 10 _1 Wm _1 K _1 . Experimen- 
tal conditions of P = 1W, R = 0.5mm give a temperature of 
AT « 100K. With dn/dT «3x 1(T 4 , n = 1.5 a deflection 
angle of « 3° is expected. When D is chosen large enough, 
one may easily observe the deflection Ax = D tan (0) on a 
screen. 



where the introduced interaction parameter v = ^ de- 
notes the strength and polarity of the potential. A pos- 
itive value of v > corresponds to a repulsive, and a 
negative v < to an attractive potential. The analytical 
solution, first given by Gordon in 192& 43 , to the equation 
is achieved by the ansatz (r) 

fY 2 



e f (r — z), wherein 



p 2 + z 2 . The complex- valued function / describes 
the perturbation of the incoming plane wave. Inserting 
this ansatz and the Laplacian in cylindrical coordinates 
into the SE leads to 



2uk 



f(r-z) = 0. 



' d 2 ld_ 2 . k d_ _d 

dp 2 p dp dz dz 2 r 

(24) 

Defining the function g (x) by / (r — z) = g (x) with the 
variable substitution x = ik (r — z) one obtains Eq. (|25|), 
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which is the hypergeometric differential equation for g (x) 

g(x) = 0. (25) 



d2 m ^ d • 
x— + {l-x) — +iv 



dx 2 



dx 



Therefore, the solution wave function for the case of 
an incident plane wave may be written dowrP* 

* c ( r ) = e -*"e ikz T (1 + iv) 1 F 1 (-w; 1; ik (r - z)) 

(26) 

The pre-factors ensure a normalization to unity \^c\ 2 — 
1 at large distances z — » oo from the scatterer. Also, 
the wave-function reduces to the incoming plane wave 
for vanishing perturbation, i.e. S&c ( r ) = e %kz for v — 0. 
An asymptotic expansion of the confluent hypergeomet- 
ric function for large kn = k (r — z) — » oo allows the wave 
function to be separated into a scattered spherical wave 
with angle- dependent amplitude / and a plane wave re- 
sembling the form e lkz +f (6) e %kr /r, although both terms 
will include logarithmic phase distortions due to the long- 
range character of the Coulomb potential 9 ! Apart from 
corrections vanishing for r — >• oo, the scattering cross- 
section reads 



da 



sin 4 (0/2) 



(27) 



On the positive z-axis the wave- function Eq. (26) 
satisfies^ |* c (z) | 2 = 2ttv / [e 27ru - lL The amplitude 
of this wave function is shown in Fig. (pi. 

Now we will write down the scalar Helmholtz-Equation 
for lightP, \7 2 U + k 2 [n (r) 2 /n%] U = 0, with the refractive 
index profile given by Eq. d5l. One should think of this 
equation as being the analogon to the SE with non-zero 
energy E = n\/2 and potential energy V = -n(r)72 + 
tvq/2 as before in Section IV We find: 



v 2 /7 + r 



l 



2AnR 

n r 



An 2 R 2 
nlr 2 



U = 0. (28) 



A comparison with Eq. (23) allows the identification of 



the interaction parameter v as 



AnR 

n 



k 



(29) 



to first order in the small quantity An /no <C 1. In- 
fact, this identification could have been guessed without 



this inspection simply by the definition of the parameter 
v ~ §§ ano ^ ^he classical correspondences found earlier 
with its prescription 2E/C — » £. The solution Eq. 
(26), to the SE of the Coulomb scattering problem may 



thus be used to find the scalar field amplitude U in the 
case of diffraction by the inhomogeneous refractive in- 
dex field Eq. (J5|. Again, the particle problem may thus 
be used to obtain results for its corresponding optical 
phenomenon, similar to the mechanical-optical analogon 
which is the "F=ma" -optics framework. The strength- 
and polarity-encoding parameter v is found to be pro- 
portional to oc which was the parameter describing 
the ray-trajectories. 

Another view on this equivalence is obtained by look- 
ing at the Kirchhoff diffraction for this refractive index 
profile. We have recently used the Fresnel- grade approxi- 
mation of the diffraction on this refractive index profile to 
describe the photothermal signal of single heated nano- 
particles^4 We may write for the scalar field amplitude 
U in the image-plane located at a distance z behind the 
aperture-plane the following diffraction integral: 

u= k_ eikz+i ^j\ ei ^ Jo (kp^ exp { . Axp) pdp 

(30) 

Here, U a is the field in the aperture-plane at z = 0. 
The collected phase advance Ax p may be computed in a 
straight-ray approximation to yield, neglecting an addi- 
tional constant phase, 

A Xp (p) = k f n (V^ 2 +P 2 ) dz « -2k R An In (g) , 

and is caused by a travelled distance I both in front and 
behind the lens. To mimic the quantum mechanical ini- 
tial state of a plane wave used before, we will consider 
a unit amplitude plane wave in the aperture plane, i.e. 
U a = 1. The abbreviations k = /co^o? a = AnR/ no 
and £ = — ^ are used from now on. We will also 
use the relations iFi (a; b\ — x) = e~ x \F\ (b — a; 6; x), 
log(— i) = —m/2. While in our previous paper we have 
studied the limit of z <C kuo 2 , we will now retain the 
^-dependence. Using the symbolic arithmetics program 
Mathematical or the integral representation of the conflu- 
ent hypergeomeric function \F\ in combination with the 
Bessel function representation through q^i, one finds the 
image-plane amplitude to equal 



J 



U Or, z) = ^e ikz + ih £ |C" 1+<fco r (1 - ika) X F X (l - ika; 1; -£fj) 



■ t kx z 

' 2z 



(31) 
(32) 



I 

Which already resembles the exact solution S&c found metric function in the forward direction using the relation 



earlier. If one considers the argument of the hypergeo- 



2 2 2 



2z [r — z], one may write i- 



ik [r — z], 
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such that an agreement in form is found between Eq. ( 32 ) 



and Eq. (26) up to a logarithmic phase-factor. Upon in- 
spection of these two equations, we may thus also arrive 
at the identification v — >• — fca, which is the same as stated 
in Eq. (29). The absolute value of the final expression 
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FIG. 5. Top: \^ c (r,z) \ 2 - 1 from Eq 
shadow lines (green) and the x (z) = z tan 



([26} and the GOA 
( o) lines (blue) 

encl osin g the "near" zone. Bottom left: Zoom. Bottom right: 
Eq. ( [32| . The first line-scan position of Fig. [6] is indicated by 
the dashed line, tan (6) = x/z, v — 0.214, A = 635nm. The 
black lines in the bottom half depict several ray trajectories, 
Eq. ( |12[ ). Upper half shows trajectories for 8x enhanced lens- 
strength, i.e. v — 1.71. 



for the field amplitude shows that for An = a = one 
obtains U(r) = e zkz , i.e. the plane wave emerges unper- 



turbed (1F1 (a; b;x) — 1 as x 



We have thus 
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for 



demonstrated the equivalence of the plane-wave (pw) 
quantum mechanical Coulomb scattering problem and 
the pw diffraction by our specific thermal lens n(r), to 
first order O (An/n). 

The connection shall now be demonstrated between 
these wave-mechanical descriptions and the previously 
studied classical cases. While the shape of the wave- 
function amplitude already resembles the family of tra- 



jectories of a given energy but varying impact parame- 
ter, i.e. Fig. [T}d, the resemblance is misleading. For a 
given constant wavelength A and thus constant wave num- 
ber, the diffraction pattern is independent on the mag- 
nitude of the thermal lens, i.e. An, while the family 
of trajectories would change. For plane wave scatter- 
ing/illumination the spatial features and patterns of the 
perturbed wave- amplitude do not resemble the trajecto- 
ries and shadows in the near field as predicted by ge- 
ometrical optics^. While classical dynamics and scat- 
tering descriptions require the notion of paths and tra- 
jectories, in the wave-mechanical scattering description 
no clear correspondence exists for the case of plane wave 
or wide beam scatterimJ^. It is only in the far field 
that the classical average particle number density co- 
incides, up to an additional zero- mean oscillation with 
an undetectably high spatial frequency , with the QM- 
wavefunction amplitude^ and thus with the classical ex- 
pressions for the cross-section (dcr/dfi). We will there- 
fore now formulate the correct limit which connects both 
wave and the classical descriptions. 



B. QM: Wave Packet Scattering 

In order to draw the connection between the classical 
and the wave pictures both in optics and quantum me- 
chanics as we have described above, it is necessary to 
reconsider what constitutes this classical limit such that 
a recovery may be demanded. In the previous paragraph 
it was already shown that the plane- wave approach does 
not resemble the classical pictures apart from the total 
far-field scattering-cross-section. In optics the transition 
is reached by letting the wavelength go to zero, A — » 0. 
Then, the wave-front normals will follow the trajectories 
described by ray-optics 15 . Due to the possession of the 
exact solution in the case of quantum mechanical wave 
theory, its transition to the classical particle trajectory 
picture will be outlined here. Here, we will investigate 
the QM scattering problem of a wave-packet (wp) and 
will afterwards instead of following phase-front normals 
quantify the scattering of a specific wave-packet which 
will resemble a confined minimally spreading ray. Fol- 
lowing Baryshevskii et alP and similar works=™^, one 
may write an initial wave packet localized near (i.e. fo- 
cused at) r at time t = as: 



*7(r,0) 



J dkA(k) e ik '( r - r °) 
A(k) = J drG(r)e- ir '( k - k °). 



(33) 
(34) 



Such a superposition may even have the familiar phase 
anomaly known from wave optics as the Gouy phased 
The functions A (k) and G (r) define the wave-packet 
form in momentum- and real-space, respectively. For dif- 
ferent momenta k and thus possibly different interaction 
parameters = or = k^~ x the solution formerly 
written down for k 



fez, Eq. (26), will now be given 
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in a fixed coordinate frame for arbitrary direction of the 
incident wave- vector: 



XZ— plane 




FIG. 8. Initial wave-packet amplitudes \^q P (r) | . Im- 
ages show two focused WPs with zo = +400nmz and 
a?? = {15°, 30°} degrees. Axial Gaussian fits yield cjo = 
23nm + 0.8cj^ with uj$ = 2/[ka#] and axial Gaussian widths 
cjo,2 = 120nm + 1.32zr with zr = kcj#/2. 



^ (r) = e-^ fc e ik r r (1 + iu k ) 1 F 1 (-iu k ; 1; i (fcr-k-r)) . 

(35) 

The time evolution of an arbitrary initial wave-packet as 
described by Eq. (37) will be determined by the superpo- 



sition \I/^ P of the individual plane-wave solutions corre- 
sponding to the pw-spectrum components of this initial 
wave-packet: 



*S P (r,t) = |dkA(k)e- 



ikr ° (r) e 



(36) 



Now, assuming that only different incident angles will 
contribute, we have a constant momentum magnitude 
of |k| = k. Further, we will consider an azimuthally 
symmetric angle distribution for the wave-vector spec- 
trum such that it takes the form A (k) = A (fc, tp) = 
5 (k — k) A (#) /27r, where the wave- vector has been ex- 
press in spherical coordinates. The initial wave-packet 
then reads: 



FIG. 7. Geometry for wave- vector k of (r) used in the 
wave-packet superposition (r, £). The azimuthal angle of 
the xz-plane is cj) — for x > and 4> = tv for x < 0, while 
the polar angle is = arccos (z/r) with r 2 



2 i 2 
X + Z . 



*7(r,0) 



d$ A (l?) e - ik izo-rcos(0)] cos(tf) ^7) 

x sin (#) J (fcr sin (0) sin (#)) 



More specifically, we will choose the polar-angle spectrum 
A (#) to be a Gaussian with an angular width of cr#, 
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A (i?) = exp 



(38) 



which results in a focused wave-packet that has similar 
properties as a TEMOO-mode Gaussian beam with a char- 
acteristic width-scale given by uj$ = 2/[ka&] (see Fig.[8|. 
This immediately implies that the angular spreading of 
the WP decreases, i.e. becomes paraxial and resembling a 
ray, when the wavelength decreases since <j# oc \/uj$. In 
fact, as is the case for the beam emerging from a laser- 
pointer, its lateral intensity follows closely a Gaussian 
distribution of width cjo, i.e. \^q P (x) | 2 oc exp (— 2x 2 /uJq) 
while the axial intensity pattern is enveloped by a Loren- 
zian profile with a corresponding of range zr = fccJo/2, 



i.e. 



ocl/ [1 



2 /z%\ (solid thick lines). A fit 



in axial direction may also yield some axial Gaussian with 
widths uoq^z which fit the central bumps (thin solid lines). 
The phase-pattern also shows the Gouy-phase anomaly, 
i.e. a phase advance of tt as compared to a spherical wave 
emanating from the focus. 



The solution to the time-dependent SE, Eq. (36), with 
the specific choice of the initial wave-packet as described 
by Eq. (38), then reads: 
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*c p M) 



■- \ dtf / d<, 

JO JO 



dip sm(tf)A(tf)e- ikzoCOsW e-% u e ik - r r(l + iv) 1 F 1 {-iv; 1; i (kr - k • r)) e 



■hk±. 
2m 



(39) 



with the dot product of the wave-vector and the ra- 
dius vector in spherical coordinates k • r = k r cos (a) = 
k r [cos (6) cos (#) + sin (0) sin (#) cos (0 — y?)] (see Fig. 
[7]), z/fc = z/ and k • r = k zq cos ($). 

To obtain the ray-limit, we finally choose beams of 
a finite width uj$ and small angular spread (i.e. parax- 
ial) with an lateral offset in x-direction of the resulting 
stretched wave-packet to some zq ^> In this case 
one must set k • r = k zq sin (1?) cos (tt — cp) (see Fig. 
[7]). We then find the wave-packets to be distorted by 
the scattering process such that its probability ampli- 
tude |\I/^ P (r) 1 2 follows the classical Rutherford scatter- 



ing trajectory r(0), Eq. (20) with the plane polar angle 



4> now being the polar angle with the impact param- 
eter set to the initial WP lateral offset b —> zq and the 
lens strength parameter £ — >• kjv. This correspondence 
is shown in Fig. [9] and is the analogue to Fig. [TJd. This is 
the expected classical property of the quantum mechan- 
ical scattering description^. 



1.0 ^ 




-1 1 

z [jLvm] 



FIG. 9. Initial wave-packet amplitudes \^q P (r) | 2 . Im- 
age shows two focused WPs with z — — 0.600/xmx and 
zo = 1.0/xmx with a spreading of a$ = 3° degrees and a cor- 
responding width-scale of oj$ — 0.13/xm. The strength of the 
potential was v — 17.1, the wave-number was k — 289/un -1 . 
The scattered WPs follow the classical photonic Rutherford 
trajectories (red-black dashed), similar to Fig. [TJd, and avoid 
the shadow region (green line, textured area). 



C. Perturbation quantification: The normalized 
difference / photothermal signal $ 

In a photothermal microscopy experiment a focused 
beam is used to probe the refractive index lens of a heated 
nano particle. A Gaussian beam provided by a TEM00- 
mode operated laser will be focused by a microscope ob- 
jective. It will thus have a total angular spreacP^of twice 
the divergence half-angle #div = ^j- For A = 635 in a 
medium with no = 1.46 and a beam waist ujq = 281 nm 
one finds 6aw ~ 28°. Therefore, one needs to consider 
both wave- mechanical pictures with focusing. First, we 
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FIG. 10. Focused wave-packet scattering on a Coulomb po- 
tential with no axial offset of the initial wave-packet and finite 
spreading width a$ = 15°. Depicted are |^^ p | 2 (top) 
(center) and their normalized difference (bottom). 



Having previously established the connection to the 
optical wave-mechanical framework, one may also think 
of these trajectories to be refracted narrow gaussian 
beams. The power of the general beam description is 
shown in the following via its application to photother- 
mal microscopy (see also Fig. [l0|). 



will discuss the QM analog on to the rel. PT signal. To 
this end we will need to focus on the narrow forward 
direction interference-zone which is otherwise neglected 
in standard treatments of Coulomb scattering problem 
as its angular extent shrinks to zero at large distances. 
However, as it will turn out, it is exactly this interference 
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which causes the PT signal and the interference-zone is 
expanded to detectable angular extents via the finite- 
width plane-wave wave- vector spectrum. In the small an- 
gle interference domain (see Fig. (J6|), we need the asymp- 
totic form of z), Eq. (26), for small r] — r — z. As 
stated before, in the treatment of QM Coulomb scatter- 
ing the opposite limit of kr] — >• oo is usually considered^. 
Now, series expansion of the confluent hypergeometric 
function reads iFi(a;b;z) = 1 + + G[z 2 ). Using 
e x « 1 + x for small x <C 1, Eq. (26) thus becomes to 
order O (rj) and to order O (v): 



(r) 



1 



ikz 



[1 



ij E is] [1 • 
\- kr] + ii 



iv x ikrj] (40) 

/c 2 7? 2 ' 
-7^; + — r- 



since 77 = r [1 — cos (0)] = 2r sin 2 (0/2) we have to first or- 
der in v again the squared modulus of the wave-function: 

\®c(0,r) | 2 « l + i/ [4/crsin 2 ((9/2) - tt] (41) 

This approximation nicely fits the central bump. For 
zero angle, i.e. in forward direction, the finite value 
\^c (0 = 0) | 2 = 1 — ttu. A more rigorous demonstration 
of this limit can be found in loc. Eq. (42) of reference^. 
The Photothermal signal analogue would thus read, for 
plane- wave illumination and small angles: 



$ = |^c| 2 - 1 ~ -kv + krvQ 2 



It is a parabola which cuts the 



(42) 



c\ 2 = 0-axis at 60 = 
±Y^^ independent of the strength of the perturbation 
£. This is the width of the central bump and defines the 
"near" -zone as in loc. Eq. (9) of reference^ derived on 
ground of different arguments. It depends on the distance 
r to the scatterer. Since analytical progress in the case 
of wave-packet scattering was here non-feasable, we will 
now concentrate on the corresponding optical scenario to 
achieve a more general expression which will generalize 
Eq. (42) to the non plane- wave case. 



If one assumes a Gaussian beam which illuminates the 
aperture plane a corresponding substitution of U a by 



Up 



exp 



UJ 2 (Z ) 



+ ikzo + i 



kp 2 



2R C (*o) 



Kg 



into Eq. (|3l|) requires the substitution of ( = 



H — 2Rc(z ) anc ^ a P re_ f ac ^or including a phase factor 
(Gouy-phase) and an amplitude factor Uq/uj (zo) 2 . The 
beam-waist of the focused Gaussian beam is denoted by 
uo (zo) 2 = Uq [l + Zq/zr] and the local radius of curva- 
ture by Rc (zo) = zo [l + z 2 R j z 2 \ and the Gouy-phase is 

Cg = arctan(z /^)- 

In the case of the diffraction formulation we can write 
for the relative photothermal signal $ on the z-axis in 
case of an illumination Gaussian beam from Eq. (31): 



$= ^=^ =e ^arg(C)| r(1 + ^ ) |2_ 1 (43) 

\Uv=P\ 



The number \v\ = kR\An\ < 1 is small, such that 
r (1 + iv) « 1 — ijev + O [y 2 ) where 7^ is Euler's con- 
stant. This means that to first order in O (y) one may 
write 



$ = 2v arg (£) = 2v arctan 



-zo 
zr 



(44) 



The plane-wave limit may be read off directly from Eq. 
(32), <I> = — 7tza It also agrees in value with the above 
Eq. (44) for large offsets of the beam-waist, i.e. \z$\ ^> zr 
and zo > (scatterer behind beam- waist). It is clear 
by the discussions and equations presented so far, that 
photo th ermal single particle microscopy as described in 
reference d 11 ^ 12 ! deals with the interference zone encoun- 
tered in Section |V A of plane- wave quantum mechani- 
cal Coulomb / Rutherford scattering. In this case the 
interference zone exhibited a vanishing angular extent 



kr 



thus being typically discard in scattering 
analysis and leading to the sin _4 -dependence of the de- 
tectable cross-section. In photothermal microscopy it is 
modified and widened up to the extend of the angular 
spread of plane- wave contributions which make up the in- 
cident beam, i.e. the angular spread of the probing laser- 
beam. Thereby, the interference zone corresponds to 
the detected angular domain of photothermal microscopy 
and determines the signal. On the other hand, if a de- 
flection is measured of the probing beam, the usual low- 
energy (wide lateral wave-packet width) limit of Ruther- 
ford scattering is attained and deflections by the scat- 
tering angle are expected in the limit of small angular 
spread and far enough lateral offsets (see discussion in 
Section [VB). 
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FIG. 11. QM wave mechanical probability amplitude (x) 



Eq. ( 26 ) at various axial coordinates z (red) and wave-optical 
scalar diffraction results, Eq. (31), for various beam waists 
(blue). The optical diffraction results have been computed 
using a Gaussian beam illuminating the aperture-plane as de- 
scribed by Eq. (43). 
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VI. APPENDIX 



A. Binet's equation for pertubative forces 



While the perturbative force here is oc 



the per- 



turbative force as obtained from the geodesic equation 30 
in the Schwarzschild metric of a massive body which en- 
ters the gravitational two-body problem is ex r -4 . Given 
below are the differential orbit equations for the per- 
turbed motion of massive particles in the potential fields 



V(r) = Cr' 1 + 
F (u) = —mL 2 z u 2 



ar 



+ fir . From Binet's equation 



d\u + u 



one readily finds 



which yields an imaginary eccentricity e s = ±i|&£|, which 
would give imaginary radii. However, the hyperbolic sine 
gives trajectories which approach from infinity if one 
sets </>o = 0. Then, the approach direction may still 
be specified by the choice of e s = — 1/ sinh (771-), how- 
ever the distance to the optical axis is then fixed to be 
b = P7 _1 tanh (^7) which does not reduce to b. The 
degree of freedom to achieve this was lost when the inte- 
gration constant 0o, s was discarded. 



D. Derivation of the differential scattering 
cross-section 



din 



2a 



3/3 



C 



(45) 



While {a ^ 0, /? = 0} corresponds to Fermat's least op- 
tical path Eq. ( pT| ), the combination {a = 0,/3^0} 
matches in form the equation of motion of a massive 
particle in the Schwarzschild metric. The null geodesic 
equation for light rays requires to put = in the 

latter one and will th us als o not correspond to the path 
described by Eq. OPEI 



Using = 2(/> — 7r and arccos (x) = arctan (Vl — x 2 /x) 
for x > we can write with Eqs. dl3|), (ll9|): 



= tt — 2|7 1 1 arctan (je) 



(49) 



Since arctan (x) + arctan (1/x) = tt/2 one has for 7 w 1, 
i.e. the first order approximation, the following relation, 
6 « 2 arctan (1/e) or equivalently 



cot [ - ] = &e. 



(50) 



B. Alternate trajectory formulation 

An equivalent trajectory formulation similar to the one 
given in reference^. The distance rc of closest approach 
for a repulsive potential for a head-on impact, i.e. b = 
0, will be used. In the mechanical case, this may be 
evaluated by setting the kinetic energy of the incoming 
particle at infinity equal to the potential energy at the 
minimum distance, E = — , yielding rc = C/E. This 
implies a photonic analog on of 2£ -1 . The trajectory 
Eqns. then read: 



^RF (0) 

b 



-^[1 
2b [ 



'(0)] 



1 



r(0) £6 



[1 + cos (0)] + sin (0) . 



(46) 
(47) 



C. The hyperbolic sine case 

Now, lets assume a solution with the hyperbolic sine 
function: 



r s (0) 



Ps 



e s sinh (7 S [4> - <f> 0iS ]) + 1' 



(48) 



again with p s = — p and 7; 



2 _ 



~7 2 . These orbits de- 



scribe trajectories which approach from infinity but fall 
exponentially fast into the coordinate origin without a 
closest distance r m . The unique infinite distance may 
be set to happen at (j) = 7r, giving the constant 0o, s = 
7r + 7~ 1 arcsinh (l/e 3 ). The eccentricity may be obtained 
here by requiring the limit of lim^^ sin ((/>) r s ((j)) = b 



The geometric definition of the differential scattering 
cross-section is 



2tt ( ^jsm(0)d0 = -27rbdb. 



(51) 



Prom Eq. ( 50 ) and the derivative cot' (x) = sin 2 (x) we 
find: 



bdb 



cot(f) . 



27 sin 2 



0\ d6> 



2 7 2^ 



(52) 



Combining Eqs. (52) and (51), and using sin (2x) 
2 sin (x) cos (x) we write: 



cot (|) 



da 



dQ sin (0) 2£ 



2/ • -2 



4e 2 sin 4 (f) 



(53) 



Only the cross section for scattering at a greater angle 
than some chosen angle, a > e, is defined, and evaluates 
to, using a further trigonometric identity for cot (|) and 

sin" 4 (|) sin (0') dfl'd^ = 47rcot 2 (f ): 



<r> e w ^ 



cot z 



e 



7T / 1 + COS (©) 



£ 2 \l - COS (9) 



(54) 



E. The GOA Shadow Region 

For the Coulomb scattering, one find d 48 * 55 -*^ 

r-RF, s [1 - cos (</>)] = 2r c , 
r s [l-cos(^)] =4C~ 1 . 



(55) 
(56) 
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In order to obtain the shadow region, one may con- 
sider for simplicity the intersection of neighboring asymp- 
totes y\ (z) and y b A db (z) and find the corresponding Z{- 
coordinate of intersection for db — » 0. There will be a 
solution for each b that plays the role of a parameter de- 
scribing the intersection-point curve yi (z) (green lines in 
the figure above). The Asymptote was Taylor-expanded 
to second order around bt; = oo to obtain approximate 
results. Using y A (z) = b + tan (6) [z + b tan (0 O — tt/2)] 



and requiring db y\{zi) one finds 



b 2 e 



Vi 0) f 



2b 



4£ 2b 2 ^ 3 



7T (2 + <) 

4(l + z£) 



(57) 
(58) 

(59) 



such that yi(z) « 2 3 / 2 ^z/l - ^ _1 7r/4. This de- 
fines a shadow-region which spans an angle # s given by 
tan(0 s ) = lim^^oo ^ (z) jz. This angle goes to zero. 



More precisely, 6 S w 2 3 / 2 



7T J_ 
4 z£ ' 
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